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Abstract 

An experimental study is conducted to determine the 
dependence of vortex generator geometry and impinging 
flow conditions on shed vortex circulation and crossplane 
peak vorticity for one type of vortex generator. The 
vortex generator is a symmetric airfoil having a NACA 
0012 cross-sectional profile. The geometry and flow 
parameters varied include angle-of-attack a, chordlength 
c, span h, and Mach number M . 

The vortex generators are mounted either in isolation 
or in a symmetric counter-rotating array configuration 
on the inside surface of a straight pipe. The turbulent 
boundary layer thickness to pipe radius ratio is S/R « 
0.17. Circulation and peak vorticity data are derived from 
crossplane velocity measurements conducted at or about 
1 chord downstream of the vortex generator trailing edge. 

Shed vortex circulation is observed to be propor- 
tional to M, a, and h/6. With these parameters held 
constant, circulation is observed to fall off in monotonic 
fashion with increasing airfoil aspect ratio AR. Shed vor- 
tex peak vorticity is also observed to be proportional to 
M, a, and h/6. Unlike circulation, however, peak vor- 
ticity is observed to increase with increasing aspect ratio, 
reaching a peak value at AR « 2.0 before falling off. 

Introduction 

Vortex generators are often used in a variety of 
fluid engineering applications where a small, but critical, 
amount of flow control is necessary for proper flow com- 
ponent performance.Common uses are in aircraft systems 
where surface-mounted vortex generators are applied in 
external (airframe related) and internal (propulsion re- 
lated) flows. Figure 1 depicts a typical application of 
surface-mounted vortex generators. In external flow sit- 
uations, such as that encountered on the wing surface in 
Figure 1, the most common arrangement is array forma- 
tions upstream of flight control surfaces where boundary 
layer attachment is often critical to flight performance. 
An array of vortex generators is also being used in the 
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aircraft inlet in Figure 1. In internal flows, vortex genera- 
tors are used to prevent excessive boundary layer growth, 
flow separation, and to reduce total pressure distortion of 
the airs tr earn ingested by the aircraft engine. These ef- 
fects occur readily in inlet ducts and diffusers due to such 
factors as duct centerline curvature, and large streamwise 
variations in duct cross-sectional area. 



Figure 1 Flow conditioning with vortex generators. 


The key to vortex generator performance is in the 
mixing and secondary velocity field created downstream 
by the shed vortex structures. If properly situated in 
the flowfield, the helical motion of the fluid in the vor- 
tex forces high energy fluid of the ffeestream or core 
flow into the slower moving fluid of the boundary layer. 
The re-energized boundary layer fluid is now more re- 
sistant to flow separation. Vortex generator induced flow 
may also be used to counter boundary layer thickening 
due to component-generated secondary flows. Recent 
experimental work in an S-duct 1 and a rectangular-to- 
semi annular diffuser 2 have demonstrated the effective- 
ness of this approach. 

The large number of parameters to consider when 
designing a vortex generator array for an aircraft com- 
ponent (such as the wing or inlet illustrated in Figure 
1) has the implication that experimental work to achieve 
optimum performance is often slow and expensive. This 
fact has motivated a few workers in computational fluid 
mechanics to assist in the optimizing process by includ- 
ing a means of representing vortex generators in their 
codes. A simple and effective way of doing this is to 
employ a model for the crossplane velocity or vorticity 
field induced by the vortex generators. This is the ap- 
proach taken in recent work by Anderson et al 3 ' 4 who 
examined multiple vortex generator array geometries in a 
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diffusing S-duct inlet using a parabolized Navier-Stokes 
(RNS) solver. A similar inclusion of embedded vortices 
in a full Navier-Stokes (FNS) code was implemented by 
Cho and Greber 5 for a constant area circular duct and a 
diffusing S-duct geometry. 

The advantage to modelling the crossplane velocity 
or vorticity distribution induced by the vortex generator is 
that a simple exponential expression accurately captures 
the crossplane velocity or vorticity distribution of the shed 
vortex. To demonstrate this, consider Figure 2. Figure 
2 illustrates the secondary velocity field of an embedded 
vortex (data on the left, model on the right) shed from an 
airfoil-type vortex generator. This vortex generator is one 
of 12 symmetrically placed vortex generators spanning 
the inside circumference of a straight pipe. The data is 
acquired in the crossplane one chord length downstream 
of the vortex generators. The model is constructed from 
a summation of 24 terms (12 embedded vortices plus 12 
image vortices) each having the form 6 : 


Ti 


1 — exp 



( 1 ) 


where I\ is the measured circulation of the ith vortex or 
image, u™ ax is the measured peak or maximum value 
of the induced streamwise vorticity field, and FU is the 
measured crossplane location of the vortex center. These 
quantities are known as “vortex descriptors” following 
the work of Russell Westphal 7 . Any such model requires 
a description of the vortex strength (here T*) and concen- 
tration (here u ™ ax , but viscous core radius, r C { , is another 
possibility). The model above is based on the classical 
“ideal viscous” or Oseen vortex. 8 

Many studies have demonstrated the remarkable 
similarities between the Oseen-based model and the cross- 
plan e structure of embedded vortices shed from vortex 
generators. Pauley and Eaton 9 demonstrated this in a 
study of delta wing vortex generators, Wendt et aL 10 in a 
study of blade-type vortex generators, Wendt and Hingst 8 


for low-profile wishbone vortex generators, and more re- 
cently, Wendt et a/. 6 on symmetric airfoil-type vortex 
generators. 

To complete the modelling of a vortex generator we 
must now show how the descriptors of the model vor- 
tex depend on downstream axial position, the geometry 
of the vortex generator, and impinging flow conditions. 
For airfoil-type vortex generators we might express this 
dependence as: 

T = r(z,c,/i, a, profile, 6, M , Re , ...), 
w mar = u> mar (z,c, h, a, profile, 6, M , Re , ...), (2) 

R = fl(z,c,/i,a,profile, 6 } M ,Re , ...), 

where “profile” refers to the shape of the airfoil, z is 
the axial position (measured from the trailing edge of the 
vortex generator), c is the chord length, h is the span, 
a is the angle-of-attack, <5 is the boundary layer thick- 
ness of the impinging stream, M is the freestream or 
core Mach number and Re is a Reynolds number. The 
vortex generator models in use at the present time only 
crudely approximate the functional dependence indicated 
in Equation 2. As a means of improving current vortex 
generator models this study was undertaken to examine 
the effect of vortex generator geometry and impinging 
flow conditions for one type of symmetric airfoil vortex 
generator (NACA 0012). This paper reports the results 
of a series of experimental tests covering a range of vor- 
tex generator h/c (or “aspect ratio”), h/6, a, and M. 
The vortex generators are mounted either in isolation on 
the inside surface of a straight pipe, or in a symmet- 
ric counter-rotating array spanning the full circumference 
of the pipe. Duct or pipe core flow conditions are sub- 
sonic with a no minal boundary layer thickness to pipe 
diameter ratio, 6/D = 0.09. Three-dimensional mean 
flow velocity data are acquired in a crossplane grid lo- 
cated approximately one chordlength downstream of the 
vortex generator model(s). This is done using a rake of 
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Figure 2 A comparison of the transverse velocity field of an embedded vortex; data on the 
left, and a model constructed from the superposition of Oseen vortices on the right 
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Figure 3 The Internal Fluid Mechanics Facility of NASA Lewis. 


seven-hole probes. A total of 14 test cases are considered. 
Comparisons of the experimental results with expressions 
derived from inviscid airfoil theory will be conducted. 

The influence of axial (z) decay was considered in a 
previous report. Reference 6. In Reference 6 the vortices 
were generated in a counter-rotating array configuration 
and axial decay was examined to a position of 17.25 
chordlengths downstream of the array trailing edge. In 
addition to the parametric study described above, the 
axial decay of a single embedded vortex is examined here 
for one condition of vortex generator geometry and flow 
variable. These results are compared to the axial decay 
results presented in Reference 6. 

Facilities and Procedures 

Test Facility 

This study was conducted in the Internal Fluid Me- 
chanics Facility (IFMF) of NASA Lewis. The IFMF is a 


subsonic facility designed to investigate a variety of duct 
flow phenomena. The facility, as it is configured for this 
test, is illustrated in Figure 3. Air is supplied from the 
surrounding test cell to a large settling chamber contain- 
ing a honeycomb flow straightener and screens. At the 
downstream end of the settling chamber the airstream is 
accelerated through a contraction section (having a cross- 
sectional area reduction of 59 to 1 ) to the test section duct 
The test section duct consists of a straight circular pipe 
of inside diameter d - 20.4 cms. After exiting the test 
section duct, the airstream enters a short conical diffuser 
and is then routed to a discharge plenum which is con- 
tinuously evacuated by central exhauster facilities. The 
Mach number range in the test section duct is between 
0.2 and 0.8 with corresponding Reynolds numbers (based 
on pipe diameter) between 0.95 and 3.80 million. Mass 
flows are between 3 and 7 kgs/sec. More information on 
the design and operation of the IFMF may be found in 
the report of Porro et aL u 



Figure 4 A cut-away sketch of the test section duct 
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Research Instrumentation and Test Parameters 

Figure 4 is a detailed sketch of the various test 
section components. 

A short section of straight pipe (labelled “inlet pipe” 
in Figure 4) connects the exit of the facility contrac- 
tion to the duct segment containing the vortex generator 
model(s). Static pressure taps located on the surface of 
the inlet pipe allow the nominal core Mach number in the 
test section to be set and monitored. 

The duct portion with the vortex generator(s) is re- 
ferred to as the “vortex generator duct”. The inside sur- 
face of the vortex generator duct (and hence the attached 
vortex generator) rotates about an axis coinciding with 
the test section centerline. 360 degrees of rotation is pos- 
sible. The rotation is driven by a motor and gear located 
in an airtight box above the vortex generator duct 

A typical vortex generator is illustrated in Figure 5. 
These vortex generators consist of an airfoil-shaped blade 
(with a NACA 0012 profile) mounted perpendicular to the 
surface of a base plug. The surface of the base plug is 
contoured to the inside radius of the vortex generator duct 
Both blade and plug are machined from an aluminum 
alloy using a wire cutting (electric discharge machining) 
process. Figure 5 includes a table listing the chord, 
span, and angle-of-attack of the vortex generator models 
examined in this study. The baseline or “reference” 
model has a chordlength co = 4.06 cms, a span ho — 1-02 
cms, and an angle-of-attack, c*o = 16°. 

To study variation in vortex generator span, chord, 
and angle-of-attack, the vortex generator models are 
mounted in isolation. To study variation in core Mach 
number it is necessary to test the vortex generators in 
symmetric counter-rotating array configurations. This 



Top view. 

Figure 5 The airfoil vortex generator. Dimensions 
are in centimeters, and angles are in degrees. 


is due to a torque-balancing problem at Mach numbers 
above 0.4; i.e. the lift force-generated torque produced by 
an airfoil mounted in isolation overwhelms the restraining 
torque provided by the turning gear and motor. Figure 
6 illustrates the array configuration used to study varia- 
tion in core Mach number. The vortex generator array 
consists of 12 blades, identical in geometry ( h = 1.02 
cms, c = 4.06 cms), but alternating in angle-of-attack, 
a = ±16°. As seen in the downstream view of Figure 
6, the vortex generators are equally spaced, circumferen- 
tially, at mid-chord, A between the mid-chord position 
of adjacent blades is 30 degrees, where 9 is the circum- 
ferential coordinate. 

The coordinate system used in this study originates 
in the vortex generator duct z — 0 coincides with the 
trailing edge of the vortex generator or vortex generator 
array. The axial location of downstream (r, 9) cross- 
planes are given in terms of blade chord, c. For example, 
the first survey location possible (with the reference vor- 
tex generator installed) is the (r, 9) crossplane located at 
z = 0.38c o . 

The duct segment downstream of the vortex gener- 
ator duct is stationary (non-rotating). This test section 
segment is referred to as the “instrumentation duct”. The 
fiowfield measurements are acquired in this duct through 
the use of a radially actuated rake-probe indicated in Fig- 
ure 4. To acquire data in an (r, 9) crossplane, the rake- 
probe is actuated over a radial segment extending from 
the duct wall to the duct centerline. The vortex gener- 
ator duct and vortex generator array are then rotated an 
increment in circumferential position, A 9, and the radial 
survey repeated In this maimer pie shaped pieces of the 
fiowfield are examined A narrow slot running the ap- 
proximate length of the instrumentation duct allows the 
rake-probe to be located at various downstream cross- 



Figure 6 The vortex generator array. 
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planes. A series of slot-sealing blocks determines the 
allowable axial location of survey crossplanes, z, (again 
with the reference vortex generator mounted) : 


Zi = 0.38c 0 , l.OOco, 10.38c O) 
14.75c 0 , 15.38c 0 , 24.75c 0 . 


The rake-probe consists of 4 seven-hole probe tips 
spaced 2.54 cms apart. These probes are calibrated in 
accordance with the procedure outlined by Zilliac. 12 The 
flow angle range covered in calibration is ±60° in both 
pitch and yaw for the probe tip closest to the wall. The 
calibration range for the outer 3 tips is approximately 
±30°. Uncertainty in flow angle measurement is ±0.7° 
in either pitch or yaw, for flow angle magnitude below 
35 degrees (pitch and yaw flow angle magnitude did not 
exceed 35 degrees in this study). The corresponding 
uncertainty in velocity magnitude is approximately ±1% 
of the core velocity, v xc . 

The circumferential extent of the crossplane survey 
grid is determined by the objectives of the study. In 
Reference 6 we were interested in studying the cross- 
plane domain of a single embedded vortex constrained 
on either side by its counter-rotating neighbors. The cir- 
cumferential extent of that domain was 30 degrees as 
determined by the geometry of the parent vortex gener- 
ator array. Surveying the full 30 degrees allowed us to 
determine quantities such as vortex angular momentum 
and transverse kinetic energy, in addition to the vortex 
core descriptors of circulation and peak vorticity. This 
data will again be useful at a later point in this paper 
when we attempt an analysis of the present results. We 
note now, however, that acquiring vortex angular momen- 
tum or transverse kinetic energy for an isolated embed- 
ded vortex would involve surveying most, if not all, of 
the pipe crossplane, an impractical requirement consider- 
ing the crossplane grid resolution employed here. Thus, 
for most of the present study, we will limit ourselves to 
obtaining only the vortex descriptors. These can be accu- 
rately measured by a survey grid covering the near-region 
of the vortex viscous core. In most instances this involves 
a grid of only 20 or 25 degrees in circumferential extent. 

Crossplane grid resolution is based both on the size 
of the vortex core and the time needed to acquire data 
with the rake-probe. In most test cases here the axial 
location of the survey grid is z = Co, the vortex core 
is approximately 1 centimeter in diameter, and is highly 
concentrated with large secondary velocities present Suf- 
ficient resolution is obtained with A0 = 1° on the grid 
interior, and Ar = 1.3 mm. The vortex core grows in 
size and becomes more diffuse as it moves away from 
the parent vortex generator. Thus for axial locations suf- 
ficiently far downstream (visited when examining the ax- 
ial decay of the single embedded vortex) we can coarsen 


the grid described above somewhat, with A0 = 1.5° and 
Ar = 1.7 mm. 


Experimental Results 


Table 1 and Figures 7 — 11 summarize the results 
of this study. Table 1 lists the test conditions for every 
crossplane illustrated in Figures 7 — 1 1, and tabulates the 
vortex descriptors determined for each crossplane of data. 
The vortex descriptors are also included in Figures 7 — 1 1. 
In Figures 7 — 1 1 the radial axis represents distance from 
the wall, in centimeters, and the circumferential axis 
represents angular position in degrees. The axial position 
of each crossplane is indicated. Note again that this is 
given in chordlengths of the model being tested and is 
measured from the trailing edge tip of the airfoil. 


The vortex descriptors originate from the transverse 
velocity data in the crossplane. The vector plots on 
the left hand side of Figures 7—11 are the measured 
transverse velocity data. The transverse velocity data is 
first converted to streamwise vorticity data following the 
relation: 

6v q vs 1 &v r 

u z - -7— H TV? 

or r r ov 


(4) 


where (u r , vo) are the transverse components of velocity 
in the radial and circumferential coordinates, respectively. 
Finite difference formulas are used to represent the spatial 
derivatives in Equation 4. The resulting streamwise vor- 
ticity fields are plotted on the right hand side in Figures 
7 — 11. Solid contour lines represent negative or “core” 
vorticity, dashed lines are positive or “secondary” vortic- 
ity. The contour increment for solid lines not labelled is 
—3000 sec" 1 . Umax is located at some grid point having 
coordinates ( r ci 0 c ). The vortex circulation T is calcu- 
lated by first isolating the region of core vorticity in the 
data field. This is done by referring to the contour plots 
of vorticity in Figures 7 — 11. A path enclosing the re- 
gion of core vorticity is defined. The outer boundary of 
the core is taken to be the location where streamwise vor- 
ticity is 1% of Umax * The circulation is then calculated 
according to: 


J V • ds, (5) 

path 


where V is the velocity vector in the crossplane, and 
s refers to the path coordinate. By using closed paths 
composed of line segments in the r or 0 coordinate di- 
rections the circulation is easily determined. Uncertainty 
estimates for all listed descriptors are given in Table 1. 
These are derived by combining the uncertainties in mea- 
sured velocities and probe placement in accordance with 
the procedure outlined by Moffat. 13 
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Table 1 - A summary of test conditions and shed vortex descriptor results. 


a) 1.00c 


9 degrees 


0 degrees 


9 degrees 



(d/2-r) 


Center of Vortex: r = 93 ems, theta = 91 deg, Circulation = -0.64 m2/sec. Peak Vortidty = -33189 (1/sec) 


b) z=235c 


9 de g rees 


9 degrees 


9 degrees 



Figure 7 Velocity and stream wise vortidty results for the axial 
decay of an isolated embedded vortex. (Continued ... ) 
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Center of Vortex: r = 9.1 cms, theta = 78 deg, Circulation = -0.46 m2/sec, Peak Vorticity = -7598 (1/sec) 


Figure 7 Velocity and streamwise vorticity results for the axial decay of an isolated embedded vortex. 
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a) aJpha=8 deg, z= 1.00c 


0 degrees 
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c) alpha=20 deg, z= 1.00c q degrees 

95 


0 degrees 


0 degrees 



(d/2-r) cms \ . ; ' (d/2-r) cms \ (d/2-r) < 

4.o\\ • / / ajo\ 4.0- 

Center of Vortex: r = 93 cms, theta = 90 deg, Circulation = -0.94 tn2/sec. Peak Vortidty = -34055 (1/sec) 


Figure 8 Velocity and stream wise vortidty results for variation in airfoil angle-of-attack. 




a) AR-2-546, Z=530c Q degrees 0 degrees 0 degrees 
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Center of Vortex: r = 93 cms, theta = 93 deg, Circulation = -0.24 m2/sec, Peak Vortidty = -15996 (1/sec) 

Figure 9 Velocity and stream wise vortidty results for variation in airfoil aspect ratio. (Continued 
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b) AR-1.910, £=4.00c 0 degr ees 0 degrees 0 degrees 




d) AR=0.849, z=150c 9 degrees 0 degrees 0 degrees 



Figure 9 Velocity and stream wise vortkity results for variation in airfoil aspect ratio. 


a) h=0.25 cms, Z=5.50c 0 degrees 0 degrees 0 degrees 



Figure 10 Velocity and stream wise vorticity results for variation in 
airfoil span-to-boundary layer thickness ratio. (Continued . . . ) 
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0 degrees 


0 degrees 


0 degrees 


b) h=031 cm s, z=2.50c 




Figure 10 Velocity and stream wise vorticity results for variation in airfoil span-to-boundary layer thickness ratio. 


a) Mach =0.25, z=1.00c 
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Center of Vortex: r = 93 cms, theta = 92 deg, Circulation = -0.72 m2/sec, Peak Vorticity = -35714 (1/sec) 
b) Mach =0.40, z= 1.00c q q degrtts 0 degrees 
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Center of Vortex: r = 9.3 cms, theta = 92 deg, Circulation = -1.15 m2/sec. Peak Vorticity = -56831 (1/sec) 


Figure 11 Velocity and stream wise vorticity results for variation in core Mach number. (Continued 


) 
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C) Mach =0.60, z= 1.00c 0 degrees 0 degrees 0 degrees 



Figure 11 Velocity and streamwise vorticity results for variation in core Mach number. 


Contours of primary velocity ratio, v z / v cz , are pro- 
vided in the middle plots of Figures 7 — 11. v cz is the 
core velocity of the pipe and is listed in Table 1. An 
additional feature of the plots in Figures 7 — 1 1 is the in- 
clusion of the crossplane profile of the vortex generator 
on the plots of streamwise vorticity (shaded region). The 
comer marked with the dashed cruciform is the trailing 
edge tip. 

Table 1 is organized into groupings of test cases cov- 
ering 5 variations in test conditions. The first four group- 
ings examine variation in z, a, aspect ratio, and k/8. 
Figure 7a, the single embedded vortex at z — 1 chord 
and shed from the reference vortex generator, is common 
to all 4 groupings and is thus referred to as the “baseline” 
test case. This is indicated in Table 1. The fifth grouping 
in Table 1, variation in core Mach number, is anchored 
by results depicted in Figure 11a (M = 0.25). The vor- 
tex generator array results in Figure 11a are borrowed 
from Reference 6. 

The general nature of the flowfield is illustrated by 
the three plots in Figure 7a. A concentrated circular 
vortex structure dominates the flow near the wall, as seen 
in the transverse velocity and streamwise vorticity plots of 
Figure 7a. The vortex interaction with the boundary layer 
is depicted in the contour plot of primary velocity ratio. 
In the downwash region of the flowfield (to the right 
of the vortex) the strong secondary flow is convecting 
a “tongue” of high streamwise momentum fluid under 
the vortex core, with subsequent boundary layer thinning. 
Conversely, in the upwash region of the vortex (on the 
left), the boundary layer fluid is being forced away from 
the wall, thereby increasing the boundary layer thickness 
here. A region of primary velocity deficit is also observed 
to coincide with the center of the vortex core. 

Decay of the Single Embedded Vortex 

Figures 7a-f illustrate the development of the single 
embedded vortex over the axial range: 1.000 < z/c < 
10.375. Figure 12a illustrates the axial decay of circu- 
lation, and Figure 12b the axial decay of peak vortic- 
ity, for both the single embedded vortex and the sym- 


metric counter-rotating array. The descriptors are non- 
dimensionalized with the values at z = 1 chord. The 
uncertainty in the descriptor measurements is indicated 
by the normal dimension of the plotted symbols (circles 
and squares). The decay of the single embedded vortex 
occurs without the influence and interaction of counter- 
rotating neighbors and so the vortex trajectory followed 
by the isolated vortex is different, leading to a different 
decay behavior. 



Figure 12 The axial decay of circulation and peak 
vorticity for both the isolated and array vortex. 


Consider, first, the development of circulation. One 
mechanism of circulation decay is wall friction, which 
sets up a spanwise component of wall shear stress oppos- 
ing the rotation of the vortex core. In the array configura- 
tion, the vortices are observed to form “upflow” pairs, as 
depicted in Figures 12b. The upflow pairs tend to convect 
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each other away from the wail thereby decreasing circula- 
tion losses through wall friction effects. In contrast, the 
isolated vortex has little tendency to move away from 
the wall (the small amount of movement observed is due 
primarily to the growth of the viscous core) and so cir- 
culation decay is greater for the isolated vortex over the 
axial range 1.000 < z/c < 10.375. 

For vortices in the symmetric array u> max reaches 
its greatest value at zjc — 1.000. We assume roll- 
up of the shed vortex is complete at this axial location 
and use z/c — 1.000 as the reference survey location 
when studying the influence of vortex generator geometry 
and flow conditions (as described below) on the vortex 
descriptors. Downstream of z/c = 1.000 the decay of 
v max is very rapid for both the single embedded vortex 
and the vortices of the symmetric array. This decay 
decreases rather suddenly at z/c « 3.5 for array vortices. 
Note the “hump” in the profile of Figure 12b. Again, 
this is most likely due to an interaction with the nearest 
neighbor vortex. Also note that while the decay of uj max 
for the isolated vortex is initially somewhat greater than 
that of an array vortex, it tapers off downstream of the 
z = 5c location. Downstream of the z = 8c location the 
peak vorticity of the isolated vortex is higher than that 
of the array vortex. 

Descriptors Versus Angle-of-Attack 

Figures 7a and Figures 8a-c illustrate the data for 
variation in angle-of-attack. Figures 13a-b plot the de- 
scriptors versus a. Circulation increases in proportion to 
a over the range of a examined (8° < a < 20°). Peak 
vorticity also increases in proportion to a, but only over 
the range 8° < a < 16°. It is interesting to note that 
airflow over a two-dimensional NACA0012 airfoil sep- 
arates at |a| « 16°. 14 Perhaps flow separation over the 
vortex generator is responsible for the flattening of the 
peak vorticity profile at a = 20°. 

Descriptors Versus Aspect Ratio and 
Span-to-Boundary Layer Thickness Ratio 

The aspect ratio of the vortex generator is var- 
ied by changing the chordlength of the model. Figure 
7a and Figures 9a-d illustrate the data for variation in 
chordlength. Aspect ratio, AR, is defined as: 



Figure 13 Circulation and peak vorticity 
plotted against airfoil angle-of-attack. 


location. We can estimate the value of the descriptors in 
this grouping (Figures 9a-d) at the 1 chord location by 
interpolating the decay behavior of the single embedded 
vortex discussed earlier. These estimates are done by 
linear interpolation and are listed in Table 1. Figures 
14a-b plot both the raw and interpolated descriptor data 
versus chordlength (aspect ratio). 

The span-to-boundary layer thickness ratio is varied 
by changing h while holding AR constant Thus when h 
is decreased, for example, c is also decreased. Figure 7a 
and Figures lOa-c illustrate the data for variation in h/6. 
Again, because the chordlength of the model is changing, 
the descriptor data has been interpolated to the z = lc 
position in Table 1. Both sets of descriptor data have 
been plotted versus h/6 in Figures 15a-b. 


AR = 4 x span/(7rc) = 4 x (2A)/(7rc) 

= 8 ft/Orc), (6) 

where span = 2 h is used to take into account the wall 
effect or “image” vortex generator. In acquiring this 
data, the probe was fixed at one axial station in the 
instrumentation duct Thus when the chordlength of the 
model changed, the axial location of the crossplane data 
grid, in terms of the model chordlength, changed as well. 
When studying the effects of model geometry variation it 
is desirable to have all descriptor data at one axial (z/c) 


Descriptors Versus Core Mach Number 

The variation in core Mach number is conducted for 
vortices shed in a symmetric counter-rotating arrange- 
ment The arrangement of the parent vortex generator 
array is illustrated in Figure 6. Figures lla-c illustrate 
the data for Mach numbers of 0.25, 0.40, and 0.60 respec- 
tively. All data is acquired at the z = 1 chord location. 
In Figure 1 la the central portion of a data grid 35 degrees 
in circumferential extent is plotted. This data grid over- 
laps the domain of the vortex captured there. If all 12 
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Figure 14 Circulation and peak vorticity 
plotted against airfoil aspect ratio. 


counter-rotating vortices shed from the array of symmet- 
rically placed models illustrated in Figure 6 are assumed 
to be of equal strength then the term “domain” refers 
to a 360°/12 = 30° circumferential sector of the pipe 
where the convective influence of one vortex is confined 
by the action of its neighbors. Measurements conducted 
over larger segments of the pipe confirm the equality of 
strength between vortices in the array. In addition to the 
vortex descriptors, integrated quantities such as vortex 
angular momentum and transverse kinetic energy may be 
calculated for the data in Figure 11a. The data grid used 
in Figures 1 lb-c cover just over 20 degrees and so do not 
survey the entire domain of the captured vortex. Figures 
16a-b plot the variation in descriptors versus core Mach 
number. Both circulation and peak vorticity are seen to 
rise in proportion to Mach number. 


Analysis and Modelling 


Circulation Descriptor 

The dependence of vortex circulation on vortex gen- 
erator geometry and impinging flow conditions can be 
approximated with an expression developed by Prandtl. 15 
For a finite span wing of elliptical planform and symmet- 
ric cross-section the circulation developed about the wing 



Figure 15 Circulation and peak vorticity plotted 
against airfoil span-to-boundary layer thickness ratio. 


in the plane containing the wing cross-sectional profile 
is: 


T = 


TTV cz (XC 

i + m 

1 + AR 


(7) 


If we assume that this circulation is turned into the stream 
(into the duct crossplane) Equation 7 becomes our esti- 
mate of the shed vortex circulation. A simple means of 
accounting for the retarding effect of the boundary layer 
is to replace v cz with v from the one-seventh power law 
profile of a turbulent boundary layer: 


--f-V 

w ’ 


and, writing Equation 7 in a more general form; 


r = 


K\Vcz&C 


G)’ 


(8) 


(9) 


Uii 
1 + AR 

The constants k\ and are determined from the circu- 
lation data using a least squares procedure: 


ki = 1.55. k*> = 0.637. 


nm 


Equation 9 is plotted against the data in Figures Ba- 
iba. Flow conditions and vortex generator geometry are 
taken from values listed in Table 1 and used in Equation 
9 in a manner consistent with the test conditions. So, 
for example, in Figure 13a, v cz is held to 81 m/sec, 
c - 4.064 cms, h = 1.016 cms, 6 = 1.778 cms, and 
AR = 0.637. We see that Equation 9 correlates the data 
well in every case. 
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Figure 16 Circulation and peak vorticity 
plotted against core Mach number. 


Peak Vorticity Descriptor 

Suppose that the mysterious behavior of peak vortic- 
ity illustrated in Figures 13b-16b is connected to a vortex 
image or wall effect A correlation for the peak vorticity 
descriptor at the z = 1 chord location is now developed 
below following these further assumptions: 

1. The characteristic moment in the (r, 6) crossplane 
can be equated to the rate of angular momentum 
production by the shed vortex. This moment is 
taken to be L • t , where L is the lift force acting on 
the vortex generator, and t is the airfoil thickness, 
t = 0.012c. The angular momentum, H Vy of the 
shed vortex about its center of rotation is determined 
from: 

= J f e x V dA, (11) 

domain 

where r c is the radial position measured from the 
center of the vortex, v is the transverse velocity 
vector, and “domain” refers to the entire crossplane 
of the vortex generator duct. 

2. A vortex model accurately represents the secondary 
velocity structure of the vortex (see Figure 2) at the 
z — 1 chord location. 



Fig. 17 A model of an isolated vortex in a circular 
duct using a superpostion of two Rankine vortices. 


Several additional simplifications are necessary to inte- 
grate Equation 1 1 and obtain a closed form solution sim- 
ilar to Equation 9. These are shown in Figure 17. Figure 
17 illustrates the vortex model and domain of integration 
used to calculate the angular momentum. The vortex 
model consists of two superimposed Rankine vortices. 
The Rankine model is used here so that the integration in 
Equation 1 1 may be evaluated in closed form. The vortex 
on the left is the modelled (shed) vortex, the vortex on 
the right is an image vortex positioned to represent the 
effect of the duct wall on the flowfield of the shed vortex. 
The Rankine model is a “patchwork” vortex model with 
the following definition: 

Tr . . , . . 

Vo - - — 2 ’ ( r ^ r °)> (viscous core) 

( 12 ) 

v$ = - — , (r > r 0 ), (inviscid outer field). 

2 AT 

The radius of the viscous core r 0 , the vortex circulation, 
and peak vorticity are related as follows: 


rl = 


TTWn 


(13) 


The domain of integration for Equation 1 1 are the rectan- 
gular regions enclosing the circular duct and vortex core. 
With these approximations the angular momentum of the 
vortex model is: 


Hv 
p • dz 


2 TC 


A(R , h) 


4T ' 


(14) 
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where: 


, 1V (2Rh-h 2 \,2R{h-2R) 

= (~irnr) [ W ~ arctan 


R-h 

h-2R 


2Rh 

'R-h 


arctan 


h - R 


+ 


flln 


f2h 2 -6Rh + 5R 2 Y 
\ 2h 2 — 2Rh + R 2 r 


(15) 

subject to the restrictions: (r 0 < h), ( h < R). Proceed- 
ing to the moment equation: 

E*'- = T“ : ir =t ' l= 


r r 

or C L v cz hct = — A(R, h) — 


4f 


(16) 

An expression for Cl is borrowed from inviscid (finite 
wing) airfoil theory: 


Cl — 


2/CiOr 

C 1 + iBi) 


(17) 


where “it” has been replaced with “tci” and 2 with 
“/c 2 ” as in Equation 9. Substitute Equations 6, 9, and 17 
into Equation 16, replace v cz with v in Equation 8, and 
substitute t = 0.012c to obtain: 


PoVezOch / h\ 7 

{4h + c)(A(R,h)-0 l hc)\6J’ (18) 

0o = 7-89, 0i = 1.51. 


Compare this expression to Equation 9 written in terms 
of h and c: 

„ 6.20 v cz ach ( h\$ no * 

(4 h + c) \6J ' 

Due to the approximate nature of the analysis only the 
functional form displayed by Equation 18 is of interest 
and so the particular values of the constants 0 O and 0i 
are discarded. To explore this functional form further two 
new values of 0\ are chosen, and 0 O is determined (for 
each value of 0\) by fitting Equation 18 to the baseline 
data point (Figure 7a). The two sets of constants are: 

0o{l) = 16422, /?i(l) = —HO, (20) 

and: 

0 O (2) = 920, 0i(2) = 10. (21) 

Equation 18, with both sets of constants, is plotted 
against the data in Figures 13b- 16b in a manner simi- 
lar to the procedure followed for Equation 9 and Fig- 
ures 13a- 16a We note first the proportionality between 
Umax, «, and V C2 displayed by Equation 18. In Figure 
13a we see this correlation represents the data well over 
the range 8° < a < 16°. This is also true for variation in 
v cz or M displayed by Figure 16b. For variation in AR 


and h/6 the two sets of constants in Equation 18 pro- 
duce correlations that are distinctly different. In Figure 
14b and 15b the solid lines are Equations 18 and 20, the 
dashed lines are Equations 18 and 21. Equations 18 and 
20 mimi c the interpolated data well in Figure 14b but do 
not represent the behavior of u mal versus h/6 in Figure 
15b. Equations 18 and 21 do a better job of represent- 
ing u max versus h/6 (particularly the interpolated data) 
in Figure 15b but largely underestimate u max in Figure 
14b. Different values of 0 O and 0i in Equation 18 do not 
alter this pattern of correlation. 


Summary 

An experimental study is conducted to examine the 
crossplane flow structure of longitudinal vortices shed 
from symmetric airfoil vortex generators. The primary 
goal of this study is to establish the dependence of vor- 
tex generator geometry and impinging flow conditions on 
the shed vortex descriptors of circulation and peak vortic- 
ity at an axial station 1 chord downstream of the vortex 
generator trailing edge. The airfoils have a NACA 0012 
cross-sectional profile and are mounted either in isolation 
or in a sy mme tric counter-rotating array on the inside sur- 
face of a straight pipe. The impinging flow conditions are 
subsonic, atmospheric air with a turbulent boundary layer 
thickness to pipe radius ratio 8/ R s 0.1 1 . Measurements 
of mean three- component velocities in the crossplane are 
used to derive the shed vortex circulation and peak vor- 
ticity. Vortex generator angle-of-attack is varied over 
the range: 8° < a < 20°, aspect ratio over the range: 

0. 64. < AR < 2.55, span-to-boundary layer thickness ra- 
tio over the range: 0.14 < h/6 < 0.57, and impinging 
Mach number over the range: 0.25 < M < 0.60. 

Vortex circulation is observed to increase in propor- 
tion to or, Af , and h/8. When these parameters are held 
constant, and airfoil aspect ratio is varied, circulation is 
seen to fall off in monotonic fashion with increasing as- 
pect ratio. 

The behavior of vortex peak vorticity is similar to 
that of circulation, with two exceptions: 

1. Peak vorticity increases in proportion to a up to 
a ss 16° only and flattens off thereafter. 

2. Peak vorticity rises in monotonic fashion with airfoil 
aspect ratio reaching a maximum value at AR « 2.0 
and falling off thereafter. 

The influence of a, M, h/8, and AR on the shed 
vortex circulation can be well approximated by an empir- 
ical relation derived from Prandtl’s formula (Equation 9). 

The influence of a and M on the shed vortex peak 
vorticity is represented by a correlation developed from 
a moment-balance in the spanwise plane of the shed vor- 
tex. The influence of AR and h/6 cannot be represented 
simultaneously with this correlation and is likely an indi- 
cation that something more than an image or wall effect 
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is responsible for the behavior of peak vorticity versus 
AR and h/S. 
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